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ABSTRACT 

We have investigated a sample of 5088 quasars from the Sloan Digital Sky Survey Second Data 
Release in order to determine how the frequency and properties of broad absorptions lines (BALs) 
depend on black hole mass, bolometric luminosity, Eddington fraction {L/LEdd), and spectral slope. 
We focus only on high-ionization BALs and find a number of significant results. While quasars 
accreting near the Eddington limit are more likely to show BALs than lower L/LEdd systems, BALs 
are present in quasars accreting at only a few percent Eddington. We find a stronger effect with 
bolometric luminosity, such that the most luminous quasars are more likely to show BALs. There 
is an additional effect, previously known, that BAL quasars are redder on average than unabsorbed 
quasars. The strongest effects involving the quasar physical properties and BAL properties are related 
to terminal outflow velocity. Maximum observed outflow velocities increase with both the bolometric 
luminosity and the blueness of the spectral slope, suggesting that the ultraviolet luminosity to a 
great extent determines the acceleration. These results support the idea of outflow acceleration via 
ultraviolet line scattering. 

Subject headings: quasars: general — quasars: absorption lines — quasars: fundamental parameters 



1. INTRODUCTION 

The nature and significance of broad absorption lines 
(hereafter, BALs) seen in some 10-20% of high luminosity 
quasars is still not apparent. The broad, blueshifted lines 
indicate that high velocity outflows are present in at least 
some quasars. The outflows are potentially important as 
they may enable accretion through carrying away angular 
momentum, and they may significantly chemically enrich 
the interstellar medium of the quasar host galaxy and the 
surrounding intergalactic medium. 

There are a number of potential explanations for ob- 
served frequency of BALs. In the interest of unified 
schemes, it is postulated that all high-luminosity quasars 
host the outflows that give rise to BALs. This postu- 
late is not without warrant, as quasars with and without 
BALs appear to show similar intrinsic observed contin- 
uum and emission-line properties (e.g., Weymann et al. 
1991; Gallagher et al. 1999). In its simplest interpreta- 
tion, the frequency of observed BALs is tied to the cov- 
ering factor of the outflow around the central black hole. 
That is, whether or not one observes the outflow in ab- 
sorption depends on the orientation of the central engine 
(e.g., Weymann et al. 1991; Goodrich 1997; Krolik & Voit 
1998). Early evidence from spectro-polarimetry of BALs 
bolstered the idea that BAL quasars were viewed edge- 
on (or nearly so) and the outflow was equatorial (e.g., 
Goodrich & Miller 1995). This explanation is difficult 
to reconcile in the face of new evidence from the studies 
of radio-loud BALs. The discovery of radio-loud BAL 
quasars (e.g., Becker et al. 2000) opens up the possibility 
of directly gauging the orientations of BALs. Mounting 
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evidence in the form of spectro-polarimetry and bright- 
ness temperatures (e.g., Brotherton et al. 2006, and ref- 
erences therein) is showing that BALs are viewed at a 
large variety of viewing angles (from 15° from the jet 
axis to nearly edge-on). 

Two plausible alternative explanations include unifi- 
cation in the time-domain or more complicated orienta- 
tion schemes. The former case purports that BALs are 
a short-duration (possibly episodic) phase in the duty 
cycle of the accreting black hole (e.g., Voit, Weymann, 
& Korista 1993; Becker et al. 2000; Gregg et al. 2000; 
Gregg, Becker, & de Vries 2006). Such a scenario also 
tends to connect BALs to even rarer objects like post- 
startburst quasars (e.g., Brotherton et al. 1999) and to 
more extreme objects ultra-luminous infrared galaxies 
(e.g., Sanders et al. 1988) in suggestive evolutionary se- 
quences. 

In the latter case, the geometry of the outflow, and 
hence the frequency with which it is intercepted produc- 
ing a BAL, is dependent on the intrinsic SED (or physical 
parameters) of the disk. For instance, in the conventional 
wisdom, the physical parameters that are most strongly 
tied to the presence of a BAL outflow is the Eddington ra- 
tio and the black hole mass (e.g., Boroson & Green 1992; 
Boroson 2002; Proga & Kallman 2004). Boroson (2002) 
and Yuan & Wills (2003) have shown that BAL quasars 
have strong Fe II and weak [O III] emission putting them 
at one extreme of the Boroson & Green (1992) Eigenvec- 
tor 1, which is thought to be driven to the accretion rate 
in Eddington units. In this scheme, BAL quasars are 
thought to be the more massive (hence, more luminous) 
analogs of narrow-line Seyfert 1 galaxies (e.g., Brandt & 
GaUagher 2000; Boroson 2002). 

After many decades of work, we now have the means 
to reliably estimate the fundamental physical properties 
(e.g., black hole mass, bolometric luminosity, Eddington 
ratio) of a quasar. Thus, we are in a position to ask 
whether (and how) the parameters of a BAL outflow de- 
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Table 1. Classification of SDSS-DR2 1.7 < 2 < 2 Quasars 





Pre-Fitting 


Post-Fitting 


Class 


Number 


Fraction 


Number 


Fraction 


LoBALs 


55 


1% 






mini-B ALs /B ALs 


562 


11.0% 


536 


10.5% 


AALs 


1898 


37.3% 


1813 


35.6% 


Unabsorbed 


2573 


50.6% 


2509 


49.3% 


Total 


5088 




4858 


95.5% 



Note. — The columns on the left indicate our subjective clas- 
sifications before applying our fitting prescriptions to measure the 
Mg II emission line region. The columns on the right indicate the 
number of quasars that had reasonable fits. Bad fits were typi- 
cally due to poor S/N, or the presence of intervening or associated 
absorption contaminating the Mg II emission line. 



pend on the physical quasar properties. With the large 
numbers of quasars available through the Sloan Digital 
Sky Survey, we can do so in a highly statistically signifi- 
cant manner, where only systematic uncertainties affect 
our results. Our approach here is to address and answer 
some of these questions about BAL quasars by determin- 
ing their fundamental physical properties and comparing 
them to normal quasars (i.e., to make differential com- 
parisons where systematics should not affect the result) 
and to look for correlations with BAL properties. We 
explain the details of our methodology in §2, our basic 
results in §3, discuss the results in §4, and summarize our 
conclusions in §5. We adopt a cosmology with Hm = 0.3, 
f2A = 0.7, Ho = 70 km s^^ Mpc^^ throughout this paper. 

2. DATA AND METHODS 
2.1. Selection and Classification 

Our sample comes from the Sloan Digital Sky Survey 

(SDSS; York ct al. 2000), Data Release 2 (DR2; Abaza- 
jian et al. 2004). Wc requested all 5088 objects classified 
as having broad emission lines, i.e., quasars, with red- 
shifts between 1.7 and 2.0. This redshift range places 
both the C IV A1549 and Mg II A2800 emission lines in 
the window of the SDSS spectra (which covers the range 
~3820- 9200 A). The C IV region permits us to identify 
BALs, while the Mg II region can be used to make virial 
mass estimates (McLure & Jarvis 2002). The BAL ab- 
sorption associated with C IV makes that line a bad 
choice for making mass estimates. When BALs are found 
associated with low-ionization species like Al III A1860, 
and Mg II, then the Mg II line is also compromised much 
of the time. This measurement is required for black hole 
mass estimation, which in turn is needed to determine 
the Eddington luminosity. In addition, low-ionization 
BAL quasars appear systematically reddened compared 
to other quasar classes (e.g., Sprayberry & Foltz 1992; 
Yamamoto & Vansevicius 1999; Becker et al. 2000; Na- 
jita et al. 2000; Brotherton et al. 2001; Hall et al. 2002; 
Richards et al. 2003; Rcichard et al. 2003a), making the 
determination of the continuum luminosity more uncer- 
tain and probably biased unless an uncertain correction is 
made. Therefore, we have excluded low-ionization BAL 
quasars from our subsequent analysis (leaving 5033 total 
objects). We make an additional cut below when we fit 
the Mg II region of the spectra. 
We subjectively classify objects (i.e., through visual in- 



spection) into three classes: (1) objects displaying clear 
signs of an outflow (mini-B ALs/B ALs); (2) objects show- 
ing no signs of intrinsic absorption (Unabsorbed); and 
(3) objects that have absorption near the C IV emis- 
sion line that may break up into discrete components 
at higher resolution (AALs^). The incidence of each of 
these classes in listed in the second and third columns 
of Table 1. In our subjective scheme, we adhere to the 
idea that an outflow seen in absorption should be rela- 
tively smooth (i.e., that the profile is unlikely to break 
up into more discrete components if observed at higher 
dispersion). If a profile appears too clumpy (i.e., con- 
sists of narrow components), then it is placed in the 
AAL class, if the absorption takes place near the emis- 
sion line. If the absorption has a narrow velocity disper- 
sion (FWHM <500knis^^), is clumpy and appears at 
a large blueshift, it is taken to result from intervening 
structures (which typically is also accompanied by nar- 
row low-ionization absorption lines). Since there is likely 
a continuum of velocity widths that arise from absorp- 
tion by outflows, there will be a cases where our classifi- 
cation is incorrect. To combat this, four of the authors 
(Ganguly, Brotherton, Gales, and Scoggins) have inde- 
pendently classified the spectra. Comparison between 
the authors leads to very few cases where there is any 
dispute, implying that our classifications are both uni- 
form and reproducible. Using only a single-epoch low- 
dispersion spectrum, we feel that this is the best that 
can be accomplished without resorting to more quanti- 
tative schemes. [Such schemes would require continuum- 
fitting procedures, such as template fitting or polynomial 
fitting, and arc beyond the scope of this investigation.] 

While our classifications are subjective, we feel they 
are more complete than the various quantitative schemes 
from the literature which are either insufficient, biased, 
or too contaminated with false-positives for our purposes. 
The BALnicity Index (BI) was defined by Wcymann 
et al. (1991) in order to be sure they classified only BAL 
quasars as BAL quasars with low signal-to-noise ratio 
and low-resolution LBQS spectra. For example, intrin- 
sic absorption systems that appear at high velocity (i.e., 
mini-BALs) can be excluded when using a BI criterion. 
Intrinsic absorption appearing within 5000 km s~ would 
also be excluded by a BI selection. Similarly, the Ab- 
sorption Index (Al) of Hall et al. (2002)^, while more 
liberal than the BI, still has some fairly arbitrary lim- 
its to help weed out blends of associated absorbers that 
may be intervening systems. Since we wish to test the 
properties of outflows (e.g., maximum velocity of absorp- 
tion, onset velocity, velocity width) as a function of other 
quasar properties (e.g., luminosity, black hole mass), it 
is important that our sample be roughly complete. We 
emphasize here that absorption by outflows comes from 
a continuum of velocity widths. While a small error rate 
in our classiflcations would not compromise our study, 

* The AAL abbreviation stands for "associated" absorption 
lines. Nominally, these are lines that have a narrow velocity- 
dispersion {< 500 km s~^) and appear near the quasar redshift 
(c|Ai;l < 5000km s~^). We do not adhere strictly to these crite- 
ria in our subjective classification scheme. C IV absorption that is 
narrow or sufBciently clumpy and appears superposed on the C IV 
emission line is deemed an AAL. 

^ See Trump et al. (2006) for a revised definition of this param- 
eter. 
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Table 2. Comparison to Trump et al. (2006) Classification Scheme 





Total 


Unabs. 


AALs 


BAL LoBAL 


This Work 


5088 


2572 


1898 


562 


55 


T06, AI> 


1206 


181 


420 


551 


54 


T06, AI< 


3882 


2391 


1478 


11 


1 


Breakdown 


of Trump et al. (2006) AI> classes 




high-ionization 


1158 


178 


400 


540 


40 


low-ionization 


48 


3 


20 


11 


14 


Hi 


835 


72 


186 


537 


40 


nHi 


321 


106 


214 


1 





H 


2 








2 





Lo 


41 


3 


16 


10 


12 


nLo 


3 





1 


1 


1 


LoF 


3 





2 





1 


nLoF 


1 





1 









Note. — TOO indicates Trump ct al. (2006). The fourth row, 
labeled "high-ionization," denotes the sum of the "Hi," "nHi," and 
"H" classifications. The fifth row, labeled "low-ionization," denotes 
the sum of the "Lo," "nLo," "LoF," and "nLoF" classifications. See 
text for a description of the sub-classes. 



inclusion of only BI>0, or AI>0 objects in a sample of 
objects with outflows would bias our results. 

We compare our subjectively-selected sample to the 
more objective BAL quasar sample from the SDSS third 
data release compiled by Trump et al. (2006). In total, we 
find that 562 of the 5088 objects in our sample appear to 
have spectroscopic evidence of high-ionization outflowing 
gas (i.e., with no accompanying low-ionization absorp- 
tion). By comparison, 1206 of the 5088 objects appear 
in the Trump et al. (2006) BAL catalog with the follow- 
ing classifications: Hi - 835, nHi - 321, H - 2, Lo - 41, 
nLo - 3, LoF - 3, nLoF - 1. [A "Hi" or "H" classification 
indicates the presence of C IV absorption, while a "Lo" 
classification indicates the additional presence of Mg II 
absorption. An "n" prefix indicates that the absorption 
width is narrow while a "F" indicates the presence of 
Fe II absorption.] That is. Trump et al. (2006) appear to 
find 1158 objects with evidence for high-ionization out- 
flows. At face value, it would seem that our subjective 
classification scheme is not more efficient at selecting out- 
flows. We present a cross- comparison of our classification 
scheme and that of Trump et al. (2006) in Table 2. In 
the first two rows of the table, we present a head-to-head 
comparison of the total number objects in our work and 
those objects that would have been flagged as BALs by 
Trump et al. (2006), as well as the breakdown with our 
classiflcations. In the next two rows (rows three and 
four), we further break down the Trump et al. (2006) 
numbers by ionization. In the remaining rows of the ta- 
ble, we present the full demographics using the Trump 
et al. (2()06) classification and their comparison with our 
classification scheme. We draw the reader's attention to 
the following comparisons: 

1. Of the 562 quasars with high-ionization outfiows 
(BALs) that we selected. Trump et al. (2006) only 
cataloged 551 objects. Of these, only 540 have 
a "Hi," "nHi," or "H" classification. The other 
11/551 objects (that are in the catalog) are given 
"Lo" or "nLo" classifications. The 11/562 objects 
that are not in the Trump et al. (2006) catalog are 
shown in Figure 1. Of these, six are cases where ei- 



ther the velocity limits of integration for Al or the 
continuous absorption criterion arc not sufficient. 
In other cases, the profile may be too shallow rel- 
ative to the signal-to-noise to allow a precise Al 
measurement. 

2. 181 of the objects classified by Trump et al. (2006) 
as a BAL appear to be intrinsically unabsorbed 
objects, typically with Mg II or Fe II absorption by 
intervening structures. Examples of these objects 
are shown in Figure 2. 

3. 40 of the objects that we classify as LoBALs are 
presented as HiBALs in the Trump et al. (2006) 
catalog. These objects appear either somewhat 
reddened or have broad absorption near the Mg II 
or Al III emission lines, which may affect our ability 
to carry out spectral fits of the Mg II region. 

4. For 400 of the objects classified by Trump et al. 
(2006) as high-ionization BAL quasars ("Hi" or 
"nHi"), it is not clear that the C IV absorption 
is necessarily a result of an outflow. In these cases, 
the absorption profiles seem as if they would break 
into more discrete components if observed at higher 
resolving power. The gas may be due to absorption 
by the host galaxy, or other structures related to 
the quasar /quasar environment. We place these 
objects in our AAL class. 

We also note that 1478 additional quasars in our sam- 
ple have AALs, that is narrow velocity-dispersion sys- 
tems appearing within 5000 km s~^ of the quasar red- 
shift. Some of these are likely to be intrinsic to the 
quasar central engine in some form. [Estimates range 
from > 20% (Wise et al. 2004) based on time-variability 
to ~ 33% based on partial coverage (Misawa et al. 2007).] 
None of these systems were found with the revised AI- 
based selection scheme of Trump et al. (2006). 

We conclude that, in this sample, the number of ob- 
jects with outfiows detectable in absorption is probably 
closer to our value of 562 (11%), and has likely been 
overestimated (through the inclusion of "false positives" ) 
in the Trump et al. (2006) catalog. While the efforts 
at quantitatively selecting outflows are undergoing re- 
vision (and have certainly progressed from the days of 
BALnicity), further revisions are required to decrease 
the number of false-positives. If our assessments are 
taken to be a truer reflection of outflow classifications, 
then we estimate that the Trump et al. (2006) catalog 
is 98% (605/617) complete toward finding HiBALs and 
LoBALs (consistent with their estimation), but suffers a 
15% (181/1206) rate of false-positives. Detailed scrutiny 
of subsamples of sources, such as presented here, should 
help in that goal. 

2.2. Spectral Fitting and Parameter Estimation 

Since one of our goals is to compare the physical prop- 
erties (e.g., continuum luminosity, spectral shape, virial 
mass, Eddington ratio) of absorbed and unabsorbed 
quasars, we carried out fits around the Mg II A2800 emis- 
sion line (over the rest-frame wavelength range 2000- 
3000 A) using the SPECFIT task (Kriss 1994). Our fits 
include a power-law continuum (with the convention 
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Fig. 1. — We show spectra of the 11 objects that we subjectively classify as HiBALs, that do not appear in the Trump et al. (2006) 
catalog. 



Fx ~ A~"), the I Zw 1 iron emission template (Vester- 
gaard fc Wilkes 2001), and two Gaussians (a broad and a 
narrow component) for the Mg II emission line. The fits 
were carried out on spectra that were reduced to rest- 
frame wavelengths calculated using the redshift reported 
by SDSS. 

Initially, we tried to fit all spectra with a single fitting 
prescription (as summarized in column 3 of Table 3) . We 
found that there was sufficient variety in the spectra that 
one prescription was not adequate, and we turned to five 
different prescriptions. These five prescriptions are also 
summarized in Table 3, which includes allowed ranges in 
the parameters, and parameters that were tied together. 
For example, in one prescription, we tied the wavelength 
shifts of the broad and narrow Mg II c;omponcnts and the 
Fe II together (i.e., no relative shift between them). In 
another prescription, all three shifts were allowed to vary 
independently. We also tried prescriptions with different 
ranges in various parameters like the magnitude of the 
shifts or the widths of the emission line components. 

We took the approach of fitting each spectrum using 
the all five prescriptions and adopting the model that 
provided the best fit (i.e., the lowest value). The bot- 
tom four rows of the table report the number of times 
each prescription was adopted for all quasars and for 
the three quasar classes (mini-BALs/BALs, AALs, Un- 
absorbed). In many cases, more than one prescription 
provides a good fit, with similar best-fit values and the 
actual difference in values is small. In all, we were 
able to achieve good fits for about 95% of the quasar 
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Table 3. Summary of Fitting Prescriptions 









Run Number 






Property Unit 


1 


2 


3 


4 


5 














Index^ 


[1,10] 


[0.1,10] 


[0.1,3] 


[0.1,3] 


[0.1,3] 


Normalization'' 


[1,10000] 


[0.1,10000] 


[0.1,10000] 


[0.1,10000] 


[0.1,10000] 


Fe template 












Scaling 


[0.001,10] 


[0.001,10] 


[0.001,10] 


[0.001,10] 


[0.001,10] 


Wavelength shift A 





M2B 


[-50,50] 





[-50,50] 


Width km s"' 


[900,10000] 


[900,12000] 


[900,12000] 


[900,12000] 


[900,12000] 


l\/lcy 1 1 riT"/~>tiri i^r^m T~inn 












Scaling 


rn 1 1 nl 


rn Qonnnl 


[U,OUUUUJ 


rn Qonnnl 


rn Qnnnnl 

[U,OUUUUJ 


V V Cl V Clditi Lrll Ollli L 


-28 221 


-28 221 


[-18,22] 


[-18,22] 


-98 1 021 


Width km s"' 


[2000,20000] 


[2000,20000] 


[2000,20000] 


[2000,20000] 


[2000,20000] 


Mg II narrow component 












Scaling 


[0.1,10] 


[0,30000] 


[0,30000] 


[0,30000] 


[0,30000] 


Wavelength shift A 


[-18,22] 


M2B 


M2B 


M2B 


M2B 


Width km 


[900,10000] 


[900,10000] 


[900,10000] 


[900,10000] 


[900,10000] 


Number of times chosen: 


711 


1687 


1693 


735 


32 


mini-BALs /BALs: 


65 


141 


233 


92 


5 


AALs: 


269 


644 


609 


271 


20 


Unabsorbed: 


377 


902 


851 


372 


7 



Note. — In all cases, the fits were carried out on the rest-frame wavelength range 2000— 3000 A. M2B: The parameter was tied to the 
same value as the corresponding parameter of the Mg II broad component. 
''The sign convention for the power-law index, a, is Fx ~ A~". 
'^Units for the power law are lO"^'' erg cm^^ s^^ A^^ 



spectra. 230 of the quasars could not be fit well due 
to poor signal-to-noise, or the presence of intervening or 
associated absorption contaminating the Mg II emission 
line. In the right two columns of Table 1, we report how 
many quasars in each classification were successfully fit 
with one of our prescriptions. Comparison of these num- 
bers with the pre-fitting numbers indicates that no single 
class was affected more than the others. 

We show in Figure 3 a comparison of fitted proper- 
ties (near UV spectral index, 3000 A luminosity, and 
IVIg II A2800 emission line FWHIVI) between three of our 
quasar classes, unabsorbed quasars, BAL quasars, and 
AAL quasars. We numerically computed the Mg II A2800 
emission line FHWIM from the sum of the two Gaussian 
components. Generally there is good agreement between 
the samples with, perhaps, the mini-BAL/BAL class 
having systematically broader Mg II A2800 emission- 
lines, flatter spectral indices, and higher 3000 A lumi- 
nosities. 

From the fitting parameters, we estimate two funda- 
mental physical properties for all the quasars, the black 
hole mass, and the Eddington ratio. For an estimation 
of the black hole mass for each quasar, we use the fol- 
lowing prescription derived by McLure & Jarvis (2002) 
which uses the 3000 A monochromatic luminosity and the 
Mg II A2800 FWHM: 



BH 



Ma 



= 3.37 



ALa(3000A)\ 
lO^^erg s~i J 



0.47 



FWHM(Mgll) 

km s^^ 



(1) 

Dietrich & Hamann (2004) note that this scaling law 
gives about a factor of five smaller mass than estimates 
based on the C IV or H/3 emission line widths (e.g., Kaspi 
et al. 2000). However, since our primary goal is to look 
at differences in quasar properties between BALs and 
non-BALs, this should not adversely affect our results. 



0.15 
0.1 
0.05 





0.1 



7) 

"^0.12 
o 

^ 0.08 
0) 0.06 
:^ 0.04 

C 0.02 



45 46 47 

log A L^(3000 A) [erg s"'] 



o 








S 



a. 




2000 4000 6000 8000 
Mg II FWHM [km s''] 

Fig. 3. — We compare the continuum fitting results between 
BAL QSOs (red or dark histogram), Unabsorbed QSOs (black or 
shaded histogram), and AAL QSOs (blue or light histogram). The 
distributions have been normalized to unit area. 
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Table 4. Outflow Absorption-Line Measurements 



Target 


LSI 

/I —In 
(km s ) 


^min 
/I —In 
(km s ) 


^max 
/I —In 
(km s ) 


SDSS J002344.36+143115.43 





6964 


27697 


SDSS J020608.64-080224.41 





32771 


39285 


SDSS J074851. 73+440303.60 





15421 


25505 


SDSS J085609.02+001357.71 





16512 


28419 


SDSS J102214.76+021428.80 





-4933 


587 


SDSS J150935.97+574300.56 





-2107 


-103 


SDSS J144403.96+565751.45 





19375 


23666 


SDSS J220900.66-001413.23 


425 


1469 


6677 


SDSS J222518.52-075918.45 





15453 


24780 


SDSS J233131. 91-001940.18 





-2090 


193 


SDSS J235312.78+143547.12 





1022 


26929 



We estimate the bolomctric luminosity identically for 
each quasar. The measured rest-frame flux at 3000 A 
from the fit to the power-law continuum is converted into 
an emitted luminosity using luminosity distances for our 
cosmology. We convert these monochromatic luminosi- 
ties into bolometric luminosities: 



(2) 



where A = 3000 A, Fx is the observed flux density at 
3000 A, Dl is the luminosity distance, and / = 5 is the 
average bolometric correction from 3000 A (see Figure 12 
from Richards et al. 2006). We acknowledge here that 
other, and perhaps more refined, estimates of the bolo- 
metric correction (e.g., Vestergaard 2004; Shang et al. 
2005). Use of these estimates would generally require an 
extrapolation of our power-law fit to another wavelength, 
which we wish to avoid. Since we arc making a differ- 
ential comparison between BAL and non-BAL quasars, 
this should not affect our results. Richards et al. (2006) 
estimate ~20% uncertainty in using a mean bolometric 
correction. This would have the effect of smearing out 
the distribution. In the absence of yet more refined ap- 
proaches to making bolometric corrections en masse, we 
feel this is best that can be done. 

We calculate the Eddington luminosity for each quasar 
given our estimate of the mass of its black hole: LecIcI = 
1.51 X 10'''^(M/Mo) erg s'^ (e.g., Krolik 1999, eq. 6.21). 
We then compute the Eddington ratio, iboi/.^Edd- In 
Table 7, we present our measurements of these quasars, 
as well as the two derived physical parameters (black hole 
mass, Eddington ratio), the fitting run that provided the 
best fit, and our classification. 

In addition to measurements of the continuum and 
emission-line properties of the quasar, we also wish to 
test how BAL properties depend on quasar properties. 
The Trump et al. (2006) catalog already presents vari- 
ous measurements for most (but not all) of our BALs. 
We adopt the BALnicity index (BI), absorption index 
(Al), and maximum velocity of absorption {Vmax) niea- 
surements from Trump et al. (2006) for the 551 BALs 
that appear in their catalog. For the remaining 11 BAL 
quasars that do not appear in the Trump et al. (2006) 
catalog, wc use; our measurements of BI and VmnK- These 
are listed in Table 4 (as well as the onset velocity of ab- 
sorption, Umin)- We did not measure Al values for the 11 
BAL quasars because the new definition of Al to refiect 
a true equivalent width requires knowledge of the (uncer- 
tain) shape of the C IV emission-line since the integration 



range starts at zero velocity. Trump et al. (2006) use 
a sophisticated template-fitting algorithm to reproduce 
the emission-line shapes, but this is beyond of the scope 
of this paper. A comparison of our BI and ?;max mea- 
surement methods to those of Trump et al. (2006) for a 
subsample of BALs indicates that the two are consistent. 

3. ANALYSIS 

3.1. Differences Between BAL and Non-BAL Quasars 

With measurements /derivations of the quasar physical 
properties, we can explore if there arc systematic dif- 
ferences between unabsorbed quasars and BAL quasars. 
The distributions of measured properties (Figure 3) do 
not show appreciable differences in either the general 
shape or in the mean values (Table 5). BAL quasars may 
be systematically more luminous (at 3000 A) and have 
redder NUV spectra ((a(BAL)) = 1.29, (a(Unabs.)) = 
1.55, with standard deviations of ~ 0.4), but only slightly 
. We note here that the absolute shape of the luminosity 
distribution is affected by the SDSS quasar selection cri- 
teria. The important observation here is that there docs 
not appear to be qualitatively significant differences be- 
tween the three quasar classes (Unabsorbed, AALs, and 
BALs); they are all within a standard deviation. 

Quantitative tests (like the Kolmogorov-Smirnov test) 
do show that the small differences in the distributions arc 
apparently significant. However, when dealing with such 
large samples, even small differences in the distributions 
can be manifest as having high statistical significance. 
We stress here that visual inspection of Fig. 3-4 show 
very clear overlaps, and quite similar shapes, between the 
three classifications in the distributions of all parameters. 

Figure 4 shows a comparison of the distribution of Ed- 
dington ratios, black hole masses, and bolometric lumi- 
nosities for the three samples of quasars (BALs, AALs, 
and Unabsorbed); mean values for those quantities are 
reported in Table 5. While the distributions imply that 
BAL quasars have systematically larger bolometric lumi- 
nosities than unabsorbed quasars, the distribution of the 
Eddington ratios for those two samples are quite similar 
extending down to very small values (iboi/^Edd ^ 0.1). 
This is a rather profound result given the expectation 
from the analysis of Palomar-Green quasars that BALs 
should have higher than normal accretion rates, since 
BAL quasars in the PG sample lie at one extreme of the 
Boroson & Green (1992) Eigenvector 1. In the conven- 
tional interpretation. Eigenvector 1 is taken to be driven 
by the Eddington ratio (e.g., Boroson 2002). 

One factor that could affect the robustness of this re- 
sult is the average bolometric correction of BAL quasars 
from 3000 A. If indeed the spectral shape of BAL quasars 
is different from unabsorbed quasars, one can question 
the validity of using the same (or similar) bolometric 
correction. If the mean bolometric correction for BAL 
quasars are, for example, higher by factor of ~ 2, then 
the mean Eddington ratio would be of order unity. How- 
ever, even if this were the case, we would still have a 
significant number of BAL quasars with relatively small 
Eddington ratios, given the broad distribution in Fig. 4. 

Regardless, at best it is not clear if the average bolo- 
mctric correction for BAL quasars should be larger than 
that of unabsorbed quasars. If the difference between 
BALs and non-BALs is inclination (see, for example. 
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Fig. 4. — We compare the distribution of Eddington ratio (top), 
bolometric luminosity (middle), and black hole mass (bottom) be- 
tween BAL QSOs (red or dark histogram), Unabsorbed QSOs 
(black or shaded histogram), and AAL QSOs (blue or light his- 
togram). The distributions have been normalized to unit area. 
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Fig. 5. — We show the fraction of quasars (relative to all quasars) 
with high-ionization outflows detected in C IV absorption as a 
function of Eddington ratio (top), bolometric luminosity (middle), 
and black hole mass (bottom). Lower limits arise from assuming 
that only the quasars in the mini-BAL/BAL sample arc BALs, 
Upper limits arise from assuming that both mini-BAL/BAL and 
AAL quasar samples should be counted as outflows. 



Brotherton et al. 2006, for a rebuttal of simple orien- 
tation schemes), then one might expect differences since 
how one views an axisymmetric object changes the spec- 
trum one observes (e.g., Krolik & Voit 1998). Optical 
spectra of BAL do tend to be flatter than that of un- 
absorbed quasars (e.g., Yamamoto & Vansevicius 1999; 
Brotherton et al. 2001; Tolea et al. 2002; Reichard et al. 
2003b, see also Figure 3). On the other hand, one might 
question the validity of the black hole mass scaling re- 
lations for BAL quasars, since they have been largely 
absent from reverberation mapping campaigns.^ Inaccu- 
rate black hole masses would, of course, affect our com- 
putations of the Eddington luminosities. 

In addition to the comparison in the mean values of 
the distributions, we also consider the BAL fraction as a 
function of various physical parameters. Figure 5 shows 
how the outflow fraction changes as a function of Edding- 
ton ratio, bolometric luminosity, and black hole mass. In 
each of the panels, we provide an upper limit and lower 
limit to the outflow fraction. The lower limit assumes 
that only quasars in our BAL sample should be counted 

® To date, only one BAL quasar, PG 1700-f 518, has had a re- 
ported black hole mass from reverberation mapping (Vestergaard 
& Peterson 2006; Kaspi et al. 2005; Peterson et al. 2004; Kaspi 
et al. 2000). 



as outflows: 



Outflow Frac. > 



BALs 



BALs -I- AALs -I- Unabsorbed 



(3) 



However, it is still possible with our subjective classifica- 
tion scheme that some outflows (detected in absorption) 
were placed in the AAL class. Thus we place an upper 
limit on the fraction of quasars with detected outflows 
by assuming that all objects classified as BAL or AAL 
should be counted: 



Outflow Frac. < 



BALs + AALs 



BALs + AALs + Unabsorbed 



(4) 



There are a number of possible origins for AALs (other 
than outflowing gas) so the incompleteness of the BAL 
sample is unlikely to be very large. Thus, the upper limit 
should be treated as a very conservative value. 

Generally, there does seem to be a slight trend of in- 
creasing outflow fraction with increasing Eddington ra- 
tio, bolometric luminosity, and black hole mass. How- 
ever, the distributions are, to within our conservative 
upper limits, consistent with being almost constant. 
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Table 5. Sample Mean for Fitted and Derived Parameters 



Quantity 




Unit 


mini-B ALs /B ALs 


AALs 


Unabsorbed 


(ALa(3000A)) 
(Mg II FWHM) 


10*5 erg 
10^ km s"! 


7.04±0.24 (5.67) 
1.29±0.02 (0.45) 
4.94±0.06 (1.48) 


5.88±0.14 (5.93) 
1.50±0.01 (0.48) 
4.84±0.03 (1.33) 


5.30±0.11 (5.48) 
1.55±0.01 (0.42) 
4.75±0.03 (1.28) 


(-^'bolZ-f'Edd) 

(Mbh) 

(-^'bol) 


10* M© 
lO'*'^ erg 


0.46±0.01 (0.33) 
6.18±0.21 (4.95) 
3.52±0.12 (2.84) 


0.41±0.01 (0.31) 
5.16±0.09 (3.80) 
2.94±0.07 (2.96) 


0.41±0.01 (0.32) 
4.75±0.07 (3.32) 
2.65±0.05 (2.74) 



Note. — The mean values for each fitted property, continuum luminosity, spectra index, and Mg II A2800 emission-line FWHM is 
reported for each quasar class. Parenthetical numbers indicate the standard deviations of the distributions. 



Table 6. Spearman Rank Correlation Tests 





BI 


AI 




{v) 


ALa(3000 A) 

"NUV 

Mg II FWHM 


0.047/0.286 
0.142/0.001 
-0.113/0.010 


-0.018/0.684 
0.001/0.974 
-0.069/0.113 


0.264/7.5 X 10-1" 
0.227/1.4 X 10-'^ 
-0.207/1.6 X lO-'' 


0.200/3.7 X 10-^ 
0.295/4.8 X 10-12 
-0.194/7.5 X 10-^ 


Mbh 


-0.080/0.068 
0.111/0.010 


-0.079/0.069 
0.044/0.313 


-0.042/0.331 
0.305/8.0 X 10-13 


-0.063/0.150 
0.263/8.6 X 10-10 




Width 


Max. Depth 


Num. Trough 




ALa(3000 A) 

onuv 

Mg II FWHM 


0.021 / 0.637 
0.058/0.186 

-0.055/0.208 


-0.282/4.8 X IQ-ll 
-0.125/0.004 

1.7 X 10--*/ 1.000 


0.142/0.001 
0.004/0.920 

-0.103/0.018 




-^'bol/^Edd 


-0.049/0.258 
0.051/0.241 


-0.160/2.3 X 10"* 
-0.148/6.0 X 10-* 


-0.010/0.817 
0.140/0.001 





Note. — For each entry, we list the Spearman rank correlation statistic followed by the probability of the null hypothesis. The statistics 
for the bolometric luminosity are the same as those for ALa(3000A). 



3.2. Correlations With BAL Properties 

In Tabic 6 wc report Spearman rank correlation statis- 
tics (and the corresponding probabilities of the null hy- 
pothesis) for correlations between the measured BAL 
properties from Trump et al. (2006, BI, AI, i)max, 
{v), etc.) and our measured/derived quasar properties 
[iA(3000A), ttNuV) i'boi/-^Edd, etc.]. In this correlation 
table, we have excluded the 11 BAL quasars that do not 
appear in the Trump et al. (2006) catalog. The BAL- 
nicity and absorption indices do not appear to be corre- 
lated with any of our measured/derived quasar proper- 
ties. This is not surprising given these quantities, which 
are essentially variants of an equivalent width, are com- 
plicated functions that depend on the column density of 
the outflow, the covering factor of the flow, dilution by 
scattered light, and the overall kinematics of the flow. 

The single quantity that appears to be correlated with 
most quasar properties is the maximum velocity of ab- 
sorption. Highly significant (at least in terms of the 
Spearman rank statistics) correlations are found between 
^'max and the Eddington ratio and the 3000 A luminos- 
ity. To a lesser extent, fmax slso appears to be correlated 
with the NUV spectral index. Extreme values of t;max are 
purported to measure the terminal velocity of the mass 
outflow. In the current paradigm of radiatively-driven 
outflows (e.g., Arav et al. 1994; Murray et al. 1995; Laor 
& Brandt 2002), the significance of these correlations is 
not surprising. Naively, the number of photons (i.e., 
luminosity), distribution of momenta (i.e., the spectral 
shape), and black hole mass should play a role in deter- 
mining the ability to drive an outflow. (Ionization of the 
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Fig. 6. — Wc show the correlation between the maximum velocity 
of absorption and the Eddington ratio. BAL quasars in the Trump 
et al. (2006) catalog are shown as filled circles. New BAL quasars 
from this work are shown as unfilled circles. The curve is a fiducial 
to show the predicted scaling relationship from Hamann (1998) and 
Misawa et al. (2007). The normalization of the curve is arbitrary. 



gas will also play role, and we return to this issue later.) 
Oddly enough, there is not a significant correlation de- 
tected with black hole mass, but there is with Eddington 
ratio which incorporates the black hole mass. 
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Fig. 7. — We show the correlation between the maximum velocity of absorption and 3000 A luminosity (left) and near ultraviolet spectral 
index (right). BAL quasars in the Trump et al. (2006) catalog are shown as (red) circles. New BAL quasars from this work are shown 
as (blue) triangles. Filled (green) squares are the Palomar-Green quasars from Laor & Brandt (2002), while unfilled (green) squares are 
LBQS BAL quasars from Gallagher et al. (2006). Unfilled (cyan) diamonds indicate averages of timax and L>,(3000 A)/aNUV fof the Trump 
et al. (2006) BALs in bins that are indicated by the ticks on the (blue) horizontal line. The error bars on the points indicate the statistical 
uncertainties in the mean values. The bins were chosen to have equal numbers of objects (136). The green curve is the Vniax-^'A(3000 A) 
fit from Laor & Brandt (2002), while the shaded region indicates the reported Icr uncertainty. 



3.2.1. Wmax VS. Eddiugton Ratio 

Misawa et al. (2007) and Hamann (1998) show that the 
terminal velocity of the outflow should scale with the Ed- 
dington ratio [wtorminai (-^^boi/^Edd)^^^]- In Figure 6, 
we plot fmax against the our derived Eddington ratio for 
our BAL quasars sample (BAL quasars in the Trump 
et al. (2006) sample are shown as filled symbols, while 
the additional 11 BAL quasars are shown with unfilled 
symbols). We show a fiducial curve with the predicted 
scaling and an arbitrary normalization. Overall, there 
does appear to be an upper envelope that seems to scale 
as predicted. An actual empirical fit to the envelope 
would probably favor a steeper (i.e., larger exponent) 
scaling. 

We note one BAL quasar that appears largely dis- 
crepant with the apparent envelope, SDSS J145408.25 + 
045053.54 (zcm = 1.98127). We derived an Eddington 
ratio of logLboi/^Edd ~ —1.25. In addition to the clear 
broad C IV absorption-line. Trump et al. (2006) de- 
tect a narrow C IV absorption-line system (wfwhm ^ 
500 km s^ ) at Zabs — 1-7654, and report a maximum 
velocity of absorption Umax ~ 23184 km s^ . However, 
in addition to being narrow, this system also shows ab- 
sorption in a plethora of low-ionization species (Al II-III, 
Mg II, Si II, Si IV, Fe II) and, arguably, should not be 
counted as intrinsic absorption. If this system is dis- 
counted, then Umax for this BAL quasar would be closer 
to 14000 km s~^, which is more consistent with the ap- 
parent envelope of the other BAL quasars. 

3.2.2. w,nax vs. Luminosity and Spectral Index 

In addition to the Eddington ratio, the maximum ve- 
locity of absorption also appears to be correlated with the 
3000 A luminosity and the NUV spectral index. Figure 7 



shows plots of Umax agaiust these latter two quantities. In 
both panels we distinguish between SDSS BAL quasars 
that appear in the Trump et al. (2006) catalog (plotted 
as red filled circles) and the 1 1 new BAL quasars that we 
added from our subjective search (plotted as blue filled 
triangles). In addition, we have plotted the BAL quasars 
from the Bright Quasar Survey as tabulated by Laor & 
Brandt (2002, green filled squares), and from the Large 
Bright Quasar Survey as tabulated by Gallagher et al. 
(2006, green open squares.). 

In the Wmax — L\ (3000 A) panel, we also reproduce the 
fit (corrected for our cosmology) to the soft X-ray weak 
PC quasars carried out by Laor & Brandt (2002). For 
the updated fit, we have used the same quasars as Laor 
& Brandt (2002). We note two differences in our ap- 
proach from the Laor & Brandt (2002) approach. First, 
Laor & Brandt (2002) fit all soft X-ray weak quasars in 
the PC catalog. Soft X-ray weak quasars are defined as 
those having < —2, where aox is the spectral in- 
dex derived from the flux at 2500 A and 2keV. BALs 
are a subset of soft X-ray weak quasars. Laor & Brandt 
(2002) noted that soft X-ray weak quasars appeared to 
define an upper envelope to the Vmax — La(3000 A) plot, 
consistent with the idea that the absorption does not 
always trace the terminal velocity of the outflow. We 
note here (and discuss further in §4) that, with the in- 
clusion of the SDSS BALs, it is clear that there is indeed 
an envelope (though BALs by themselves do not define 
it), and our fit should reflect that. Consequently, we 
have excluded PG 2112-1-059 from the fit since it clearly 
does not trace the envelope of the plot (see Figure 7, left 
panel where it is plotted at AiA(3000A)- lO^^^ ^erg s'^ 
and Umax ^24000 km s~^). There a number of BALs at 
smaller luminosities that have larger Umax values. Fur- 
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thermore, there are LBQS BALs at comparable luminosi- 
ties that show larger Wmax values. In carrying out the new 
fit of the envelope, we have assumed a 10% uncertainty 
in both velocity and luminosity. The parameters of the 
revised fit are: 

v = v,iL/LX^ (5) 

logLo(erg s-^) =45.0, 

a = 0.662 ± 0.004, 

logUo(km s"^) = 3.96 ± 0.29. 

We have chosen to exclude the SDSS BAL quasars from 
the fit because we do not have a convenient criterion for 
selecting which quasars define the upper envelope. Laor 
& Brandt (2002) fit all soft X-ray weak quasars (defined 
at those having < —2), since these appeared to define 
their envelope. This criterion would not work here for a 
few reasons: 

1. Since our quasars lie in the redshift range 1.7-2.0, 

most of our quasars arc not detected in, for exam- 
ple, the i?OS'^r All-Sky Survey. So we do not have 
a convenient way of computing Oqx en masse for all 
of our BAL quasars. 

2. Strateva et al. (2005) and StefFen et al. (2006) 
have quantified the dependence of on luminos- 
ity [aox = -0.1361ogL^(2500A) + 2.630]. Con- 
sequently, it is possible that soft X-ray weakness 
as defined by o^ox is not a reasonable criterion for 
defining the envelope. From the equation, the 
range in our luminosities ('~2dex) implies a scat- 
ter in the intrinsic (i.e., unabsorbed) value of cKox 
of (Ta„^ ~ 0.27. 

3. Since all of our objects are pre-selected to have 
BAL troughs (even though the maximum observed 
velocities may be small), is it possible that all ob- 
jects may be absorbed in the soft X-ray, and there- 
fore would be included in the fit if a cJox criterion 
were used. This depends on the detailed relation- 
ship between soft X-ray absorption and UV BALs. 
See, for example, Gallagher et al. (2006) for a dis- 
cussion of this. 

There is a lot of scatter in each of the panels. So, to 
bring out the highly significant correlation, we have taken 
averages of all three quantities in four bins. The bins are 
chosen to have the same number of points (136) and so 
they are equally statistically significant and independent. 
The error bars indicate the statistical uncertainties in the 
mean values. 

For the most part, if we interpret our revised Umax — 
La (3000 A) fit to the Laor & Brandt (2002) data as an 
upper envelope, the SDSS BAL quasars are very consis- 
tent with it, and fill in the entire range of velocities over 
the sampled luminosity range (which overlaps both the 
BQS and LBQS samples). Of the 536 BALs, only 2 BAL 
lies more than la above the curve, and 11 BALs arc more 
deviant than PG 1700+518, or PG 1001+054 which were 
used in reproducing the Laor & Brandt (2002) fit. 

Another BAL parameter that appears to be correlated 
with one of the quasar physical parameters is the maxi- 
mum depth of the profile (with both monochromatic and 
bolometric luminosities). We will examine the implica- 
tions of this and correlations between BAL properties 



(which is important for interpreting this apparent corre- 
lation) in a future work. 

4. DISCUSSION 
4.1. Are BALs Super- Accretors? 

From the comparison of Eddington ratios between 
BALs and unabsorbed quasars, we are led to con- 
clude that BAL quasars are not super-accretors. While 
some fraction of BALs appear to accrete near or above 
the Eddington limit (51/536 with Lboi/^Edd > 0.8), 
there arc a significant number of BAL quasars that ac- 
crete at markedly sub-Eddington rates (358/536 with 
Lboi/LEdd < 0.5). Lboi/LEdd ~ 0.5 may not be consid- 
ered a particularly "low" rate, especially in comparison 
with lower-luminosity AGN. However, it is important to 
note that BALs do not segregate themselves from other 
AGN in their distribution of Eddington ratios (Fig. 4) 
and the fraction of quasars with BALs remains signifi- 
cant down to "low" Eddington ratios (~ 0.03, Fig. 5). 

While Fig. 5 apparently shows a BAL fraction of ~0 at 
lower Eddington ratios, this is only an artifact of small 
number statistics. We only have six quasars in our sam- 
ple with Lboi/LEdd < 0.03, and one is a BAL, giving a 
BAL fraction of 16.7%. Moreover, very low luminosity 
quasars may not be capable of driving an outflow with 
a large velocity dispersion, but may still have significant 
outflows. We note that four of the six Lboi/LEdd < 0.03 
quasars do show AALs. 

This is somewhat puzzling given that supposedly BAL 
quasars lie at one extreme of the Boroson & Green (1992) 
Eigenvectors 1 and 2 (Boroson & Green 1992: Yuan & 
Wills 2003). Boroson (2002) interpreted the Eddington 
ratio and the absolute mass accretion rate as the princi- 
pal drivers of these two eigenvectors (see their Figure 7 
which provides an interpretive diagram). BAL quasars 
arc thought to occupy the high Lboi/LEdd-high A^acc re- 
gion of this diagram. Our result apparently contradicts 
this idea. 

There are two issues here in understanding why the 
Boroson (2002) interpretation may be too simple: (1) 
the BALs in the Boroson sample may not sample the full 
range of Eigenvector 1 properties allowed by the BAL 
parent population, and (2) the interpretation of Edding- 
ton ratio as the principal driver of Eigenvector 1 may not 
be correct (or complete). 

There arc only ^ 4 BALs in the sample of 162 ob- 
jects employed in the Boroson analysis. Thus, it is pos- 
sible that these four objects simply do not sample the 
full range of Eigenvectors 1 and 2 that is allowed by the 
parent population of BAL quasars. (In comparison, we 
have increased the sample by more than a couple orders 
of magnitude.) On the other hand, from an analysis of 
11 ^ ~ 2 BAL quasars. Yuan & Wills (2003) find that 
they do seem to lie on one extreme of the Fe II - [O III] 
distribution. Since this correlation is one of primary con- 
stituents of Eigenvector 1, they conclude that z ~ 2 BAL 
quasars lie at one extreme of Eigenvector 1 like their 
low-z kin. We note, however, that another ingredient 
of Eigenvector 1 is the H/3 FWHM. From Figure 2 of 
Yuan & Wills (2003), it appears that the range of H/3 
FWHM for BAL quasars overlaps completely with that 
of non-BAL quasars, so this assertion is not strictly true. 
BAL quasars at higher redshift are not simply extreme 



Physical Properties of Quasars and Outflows 



11 



Eigenvector 1 objects. 

We note that all objects in the Boroson (2002) sam- 
ple lie at low redshift {z < 0.8), while all of our objects 
lie at higher redshift (1.7 < ^; < 2). It is possible that 
there is some form of evolution in the population of BAL 
quasars (e.g., luminosity changes or differences in black 
hole mass) such that low redshift objects have systemat- 
ically high Eddington ratios while higher redshift BAL 
quasars do not. It is also possible the first principal com- 
ponent of quasar optical properties changes with redshift. 
Our result is further supported by Yuan & Wills (2003), 
who find z ^ 2 BAL quasars with Eddington ratios as 
low as ~ 0.17. This implies that accretion rate is not the 
principal driver of Eigenvector 1 properties (or that it is 
a secondary effect that is true only at low-z). 

Yuan & Wills (2003) speculate that the principal driver 
of Eigenvector 1 may instead be the availability of "cold" 
gas fueling the accretion disk as implied, perhaps, by the 
strong Fe II and weaker [O III] emission. In this interpre- 
tation, the high correlation between Eigenvector 1 and 
the Eddington ratio in low redshift PG sample is seen as 
a secondary correlation. At low redshift, most objects 
have a low fuel supply, hence low Eddington ratio, and 
the few that have an abundance of cold fuel are able to 
accrete near the Eddington rate. These low redshift, high 
Eddington ratio black holes, are those that can exhibit 
broad absorption lines. At high redshift, however, most 
quasars (i.e., both BAL and non-BAL quasars) have an 
ample supply of cold gas with BAL quasars having the 
largest supply (perhaps because they are the youngest 
sources). Consequently, all high redshift quasars are able 
to accrete near their Eddington limits, hence having sim- 
ilar (and "high") Eddington ratios. Those that have the 
largest supply of cold gas, and hence largest fueling rates, 
have the largest outflow rate and exhibit broad absorp- 
tion lines. 

We note that an important component of the Yuan 
& Wills (2003) interpretation is the differences in black 
hole mass between their high redshift sample and the 
low redshift PG sample. The high redshift quasars from 
Yuan & Wills (2003) have masses > 10^ Mq, whereas the 
low redshift PG quasars have lower masses (< 10^ M0). 
It is this fact that allows the the high redshift quasars to 
accommodate higher absolute accretion rates (since they 
have larger Eddington accretion rates) and have higher 
luminosities than the low redshift PG quasars. 

Our sample of quasars from SDSS is at similar redshifts 
to those in the Yuan & Wills (2003) sample (z 2), but 
have smaller black hole masses, smaller bolometric lu- 
minosities and span a wider range of Eddington ratios. 
Compared to the low redshift PG sample, however, we 
have larger black hole masses, larger bolometric luminosi- 
ties, but comparable range of Eddington ratios. Since we 
do not have rest-frame optical spectra for these quasar, 
their Eigenvector 1 properties are not known directly. 
Further study (e.g., infrared spectroscopy) is needed to 
understand how the Eigenvector 1 properties of these ob- 
jects relate to the low redshift PG sample and the high 
redshift Yuan & Wills (2003) sample. 

Nevertheless, if we take at face value the fact that we 
see BAL quasars in 16.7% (and outflows in as many as 
83%) of iboi/-^Edd < 0.03 objects, then we must question 
the completeness of the Yuan & Wills (2003) interpreta- 
tion. High redshift BAL and non-BAL quasars do have 



similar Eddington ratios, but not because both popula- 
tions arc accreting near their Eddington limits. Even if 
we accept that BAL quasars do indeed lie at one extreme 
of the [O III] A5007 equivalent width - Fe II/H/3 intensity 
ratio anti-correlation, we must understand first the de- 
tailed physics governing that anti-correlation and drivers 
behind it. Yuan & Wills (2003) have shown that it is not 
simply the Eddington ratio. The addition of this work 
questions if it is simply the fueling rate. There are other 
effects, such as the covering factor of the outflow and 
underlying cause(s) of that covering factor, that must be 
considered in this interpretation. 

4.2. Radiatively- driven winds? 

Our second result is a confirmation of the luminosity- 
dependent envelope to the maximum velocity of absorp- 
tion (Laor & Brandt 2002). From Figure 7, the maxi- 
mum velocities of our BAL quasars appear to obey the 
V ~ XjOm relation re-derived with the corrected 

cosmology. To explain the slope, those authors consider 
two formulations of terminal velocity and its dependence 
on luminosity: Wmax oc y^TL/R, where L is the lumi- 
nosity, R is the launching radius of the wind, and F is 
the force-multipher (Castor et al. 1975). Aside from the 

explicit \/L factor, the luminosity also comes into play 
in determining both the force-multiplier and the launch- 
ing radius. One can assume that the launching radius 
is independent of luminosity or that there is a scaling 
(e.g., R ~ VL). This gives a range of velocity scal- 

ings: fmax ~ -\/ri°'^^~°'^- The original Laor & Brandt 
(2002) best fit had a steeper slope which led the au- 
thors to conclude that there was a luminosity dependence 
to the force- multiplier. We agree with that conclusion 
with a cosmology-corrected revision of the required de- 
pendence (combining the theoretical dependence of Umax 
on F and luminosity with our revised best-fit to the 
t^max - La (3000 A) envelope): F ~ LO.32-0.82 

The force multiplier, F, relates the force from a line- 
driven (or edge-driven) wind to one that is driven purely 
by electron (Thompson) scattering. UV transitions (e.g., 
in the range 200-3200 A) are thought to be principal con- 
tributors to line-driving. The precise relationship be- 
tween the force multiplier and the luminosity (or, alter- 
natively, ionization parameter) can be complicated (e.g., 
Arav et al. 1994), and is beyond the scope of this paper. 
It is sufficient here to say that a relationship between the 
force-multiplier and the luminosity is expected (i.e., the 
two quantities are not a priori expected to be indepen- 
dent). 

Sulentic et al. (2006) employ the semi-empirical T{U) 
relationship derived by Arav & Li (1994, logF « 2.551 — 
0.536 log J7) to explicitly characterize the dependence of 

Umax with ionization parameter (fmax ^ U^-^^^). We 
clarify here that this relationship is only valid under cer- 
tain conditions which are not necessarily valid through- 
out a BAL flow: (1) the relationship does not account 
for optical depth effects'' (all lines are assumed to be 
highly optically thin), (2) the force multiplier only takes 
into account line-driving (driving from edges may also 

^ The full relation given in Arav & Li (1994) does include terms 
accounting for optical depth, but these were not included in the 
Sulentic et al. (2006) analysis. 
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be important), (3) the driving spectrum is that derived 
by Mathews & Ferland (1987) which may not be gen- 
erally applicable, and (4) the range of valid ionization 
parameters is small (—1.5 < logf/ < +0.5). A re- 
vision of this relation given in Arav, Li, & Begelman 
(1994) expanded the range of ionization parameter to 
—3 < log t/ < 1, but with a more complicated semi- 
empirical fit: logT = 3.642 - 0.1445[logC/ + 3]^ [Note 
that this still uses the Mathews & Ferland (1987) spec- 
trum and also does not account for optical depth ef- 
fects.] Furthermore, Proga ct al. (1998) and Proga & 
Kallman (2004) note that the force-multiplier is very 
sensitive to both the Eddington ratio and to the black 
hole mass. A line-driven wind cannot be formed unless 

Lho\/LEdd ^ r~^- 

Another complication in the evaluation of the force- 
multiplier (as wc have noted) and its effect on the termi- 
nal velocity of the wind is the dependence on the shape 
of the driving spectrum. From our fits, we have a hint 
of this importance in the apparent correlation between 
Vmax and Q!NUV- The correlation (Fig. 7, right panel) is 
weak, but statistically significant (Table 6). The only 
resonant lines in the NUV region of the spectrum are 
Mg nAA2796,2803 and the plethora of Fe II hnes in 
the range 2230-2600 A. In the high-ionization BALs in 
this analysis, it is unlikely the these transitions play an 
important role in line-driving. A more likely cause for 
the correlation is the loose dependence between q:nuv 
and the UV flux at shorter wavelengths, in the sense 
that "harder" values of aNuv will yield relatively more 
UV/FUV photons than "softer" values. [Note that, with 
our sign convention, F\ ^ A~", a larger value of q:nuv 
implies a harder spectrum.] We note here that, in our 
sample, aNuv is anti-correlated with La(3000A) (i.e., 
more luminous sources are "softer"). Thus the corre- 
lation between Wmax and (3000 A) contributes to the 
weakness in the corrc;lation between Wmax and aNUV- 
Consequently, the dependence of Vmax on aNUV is prob- 
ably stronger than indicated. 

Finally, we note that Wmax - ALa(3000A) curve truly 
constitutes an upper envelope. A significant fraction of 
points lie well below the curve. The details in the com- 
putation of the force-multiplier likely can explain some of 
this. However, there are a few other possibilities to note. 
As with all objects that do not show spherical symmetry, 
the orientation relative to the observer's sight-line prob- 
ably plays a role. The magnitude of this effect is not 
clear, and is affected by the fact that we do not yet have 
a full theoretical exploration of the accretion disk-wind 
paradigm. [For example, the orientation of radio-jets and 
polarization of light in the troughs in some radio-loud 
BALs seems to indicate a range of accretion-disk orien- 
tations (Zhou et al. 2006; Brotherton et al. 2006), but 
the opening angle of the wind may also be important.] 
Moreover, while the existence of the Wmax - ALa(3000A) 
envelope indicates the importance of radiation-pressure 
in driving the BAL outflow, this is not the only possi- 
ble mechanism. Magnetically-driven (e.g., Everett 2005), 
and thermally-driven (e.g., Krolik & Kriss 2001) winds 
may explain some fraction of BALs with smaller observed 

values of Wmax- 

For the Laor & Brandt (2002) points that appear at 
low Vmax, oue must also question the location of the 



absorbing gas. It is possible that some of these may 
be imrelated to the immediate quasar environment, and 
may arise from the ISM of the host galaxy, or perhaps 
other nearby galaxies. However, this cannot explain the 
objects from our sample that have small fmax- As we 
noted before, all of our objects are selected because they 
show broad wind- like profiles. Absorption from the host 
galaxy or nearby structures would have velocity disper- 
sions of a few hundred km s~^ appearing close to the 
quasar redshift, and thus would have been relegated to 
our AAL class. The low-velocity objects in our sam- 
ple may indicate the importance of non-radiative pro- 
cesses like thermal gradients (e.g., Krolik & Kriss 2001) 
or magneto-centrifugal rotation (e.g., Everett 2005) in 
driving outflows. 

A further complication is the apparently rcdshifted ab- 
sorption that we have in some of our BALs. In the case 
of host galaxy ISM absorption, this would not be a prob- 
lem since the C IV A1549 emission line is known to be 
blueshifted with respect to the systemic redshift (e.g., 
Gaskell 1983; Espey 1993; Richards et al. 2002; Vanden 
Berk et al. 2001). However, in the accretion-disk/ wind 
paradigm, the C IV broad emission and broad absorption 
lines are thought to be produced by different parts of the 
same flow (Murray et al. 1995; Baldwin et al. 1996; Mur- 
ray & Chiang 1997; Elvis 2000). [Some fraction of the 
C IV A1549 emission line may arise from different regions, 
like a virialized component (e.g., Wills et al. 1993; Broth- 
erton et al. 1994), but seems to be dominated by a wind 
component.] A full interpretation of how (or whether) 
this situation (apparently redshifted broad absorption 
sitting on top of a blueshifted broad emission line) is 
complicated by the uncertainty in what part or parts of 
the wind are being sampled by the observed absorption 
and in what determines the velocity where the emissivity 
peaks. Both of these questions rely on the orientations 
of the observer relative to the axis of the disk and to the 
direction (and opening angle) of the wind, as well as the 
dynamics and ionization structure of the wind. 

5. SUMMARY 

1. We have subjectively scrutinized a sample of 5088 
1.7 < z < 2 quasars from the Second Data Re- 
lease of the Sloan Digital Sky Survey and placed 
the objects into three classes: 562 objects show 
clear signs of an outflow; 2573 show no signs of an 
outflow (in absorption); 1898 show evidence of "as- 
sociated" narrow absorption which could have sev- 
eral different locations. The frequency of observed 
outflows is different from the recent Al-based BAL 
catalog of Trump et al. (2006) and we point out 
several differences. 

2. We have estimated black hole masses and Edding- 
ton ratios for all quasars in this sample. We find 
that there is no appreciable difference in either the 
black hole mass or Eddington ratio distribution 
between BALs and non-BALs. This implies that 
BALs in this redshift range are not super-accretors. 
Like Yuan & Wills (2003), we speculate that Ed- 
dington ratio may not be the principal driver of the 
Boroson & Green (1992) Eigenvector 1, but rather 
availability of cold gas, and the "down-sizing" evo- 
lution of black-hole mass in accreting systems. 
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3. Wc find that the maximum velocity of absorption 
as a function of luminosity has an upper envelope 
that is consistent with the best fit from Laor & 
Brandt (2002). This upper envelope is easily inter- 
preted as the terminal velocity of radiatively-driven 
wind. We find that it is also correlated with NUV 
spectral index, which may indicate the importance 
of the SED shape in governing the dynamics of the 
outflow. However, many of our BALs terminate at 
small velocities. This may indicate the importance 
of wind-orientation, or non-radiative processes in 
driving outflows. 
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Table 7. Measured and Derived Quasar Parameters 
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5.569 


3665 


2.994 


0.6159 


1 


Reg 


SDSS 


J001053.58 + 


000642.84 


51795 


0389 


348 


1 


Q'70 Q '7A 

.0(00 /U 


2.6721 


1.169 


5.664 


5655 


7.185 


0.2610 


3 


Reg 


SDSS 


J001054.41 - 


085438.66 


52141 


0651 


569 


i 




9.3966 


1.043 


16.003 


4485 


7.364 


0.7196 


4 


Reg 


SDSS 


J001111.14 + 


151034.83 


52251 


0751 


632 


i 




2.7532 


1.609 


4.573 


5700 


6.602 


A OO A A 

U.zzy4 


1 


Reg 


SDSS 


J001142.79 + 


152216. .30 


52251 


0752 


438 


i 


QQi^ 1 OA 


3.6573 


0.219 


7.918 


3270 


2.812 


r\ o o o 

0.9323 


4 


Reg 


SDSS 


J001157.76 - 


103735.21 


52138 


0652 


264 


i 


,55 iDoU 


4.2954 


0.774 


8.870 


4060 


4.573 


0.6423 


4 


A A T 

AAL 


SDSS 


J001244.62 + 


141112.97 


52251 


0752 


217 


1 
i 


, i UOOOU 


3.2471 


0.842 


5.146 


5660 


6.881 


0.2476 


4 


A A T 

AAL 


SDSS 


,10013.30.58 + 


001814.28 


51795 


0389 


385 


i 


QA /I O Q A 

,o44ZoU 


1.8880 


1.713 


.3.812 


4270 


3.401 


0.3711 


2 


Reg 


SDSS 


J001341.74 + 


143531.31 


52251 


0752 


194 


1 


.yooZou 


1.8982 


1.433 


4.437 


6480 


8.412 


0.1746 


3 


A A T 

AAL 


bUbb 


J001400.45 + 004255.49 


51795 


0389 


412 


±, 


. ( uy / ZU 


4.5593 


1.438 


7.280 


4045 


4.137 


0.5827 


3 


Reg 


bJJbb 


J001408.23 - 


085242.19 


52138 


0652 


363 


±, 


. ( 440DU 


4.5117 


0.137 


7.668 


6825 


12.068 


0.2104 


2 


BAL 


bJJbb 


J001411.00 - 


084429.15 


52138 


0652 


375 


1, 


. ( 1 zzu 


O AO AA 

2.9899 


0.287 


5.289 


5470 


6.510 


0.2690 


2 


A A T 

AAL 


bJJbb 


J001438.28 - 


010750.12 


51795 


0389 


211 


1, 


Q 1 cfi A A 


4.8011 


0.553 


9.234 


3295 


3.070 


0.9961 


2 


BAL 


SDSS 


J001453.37 - 


002827.52 


51795 


0389 


239 


1, 


.yZZDzU 


3.7010 


0.566 


8.503 


4405 


5.278 


0.5335 


2 


Reg 


bDbb 


J001507.00 - 


000800.81 


51795 


0389 


196 


1 


TAQQAA 


10.1183 


0.906 


15.970 


3020 


3.336 


1.5852 


1 


Reg 


SDSS 


J001545.33 - 


095754.12 


52138 


0652 


479 


i 


'TAC^AI A 

, /y5yiu 


1.3341 


1.949 


2.481 


4780 


3.483 


0.2358 


2 


Reg 


SDSS 


J001657.00 + 


005532.06 


51900 


0390 


331 


i 


7C;i^1 OA 


O AA/^ A 

3.0960 


1.595 


5.370 


3775 


3.123 


0.5694 


3 


Reg 


SDSS 


J001710.85 + 


135556.52 


52251 


0752 


044 


1 


O AO 1 '7 A 

,5Uoi ( U 


10.4602 


1.386 


19.863 


4360 


7.704 


0.8538 


1 


AAL 


SDSS 


J001741.85 - 


105613.30 


52138 


0652 


041 


i 


Q Afi AQ A 


3.6452 


1.407 


6.897 


3555 


3.115 


0.7331 


2 


Reg 


SDSS 


J001800.21 + 


004602.82 


51795 


0389 


611 


1 


OAAQ O A 


3.3131 


1.584 


7.342 


4690 


5.583 


0.4354 


2 


Reg 


SDSS 


J001826.80 - 


091038.78 


52138 


0652 


578 


1 


,857870 


2.9839 


0.911 


6.164 


4105 


3.940 


0.5180 


2 


Reg 


SDSS 


J001913.57 - 


103848.25 


52138 


0652 


070 


1 


,862490 


2.2511 


0.822 


4.686 


4245 


3.704 


0.4189 


4 


Reg 


SDSS 


J001959.48 - 


090809.46 


52145 


0653 


324 


1 


,792750 


1.4842 


1.450 


2.745 


6195 


6.135 


0.1481 


3 


Reg 


SDSS 


J002000.50 + 


155110.02 


52233 


0753 


374 


1, 


,751000 


1.9210 


1.737 


3.302 


4470 


3.484 


0.3139 


3 


AAL 


SDSS 


J002028.35 - 


002915.00 


51900 


0390 


234 


1, 


,926880 


2.7439 


1.765 


6.348 


4165 


4.112 


0.5111 


2 


AAL 


SDSS 


J002028.97 + 


153435.91 


52233 


0753 


430 


1, 


.764390 


4.2464 


1.156 


7.474 


4035 


4.168 


0.5938 


2 


AAL 


SDSS 


J002127.88 + 


010420.39 


51900 


0390 


443 


1, 


.819400 


4.9046 


1.518 


9.494 


7615 


16.610 


0.1893 


2 


BAL 


SDSS 


J002132.63 - 


092424.31 


52145 


0653 


422 


1, 


.869800 


2.8025 


1.081 


5.905 


5690 


7.419 


0.2636 


3 


Reg 


SDSS 


J002143.30 + 


010840.26 


51900 


0390 


534 


1, 


.901220 


2.2417 


1.618 


4.974 


4790 


4.850 


0.3396 


3 


Reg 


SDSS 


J002302.56 + 153446.39 


52233 


0753 


462 


1, 


.723460 


3.2186 


1.294 


5.268 


4790 


4.983 


0.3501 


2 


BAL 


SDSS 


J002344.36 + 143115.43 


52233 


0753 


165 


1, 


.866000 


5.4223 


1.334 


11.353 


7715 


18.544 


0.2027 


1 


BAL 


SDSS 


J002411.66 - 


004348.06 


51782 


0391 


290 


1, 


.794470 


6.6288 


1.839 


12.295 


4575 


6.770 


0.6014 


3 


AAL 


SDSS 


,1002139.61 - 


091600,58 


52145 


0653 


538 


1 


,792960 


3.5242 


1.6.38 


6.519 


1625 


5,135 


0,4204 


3 


Reg 


SDSS 


,1002118.07 - 


090055.05 


52115 


()()53 


530 


1 


,85S(i.jO 


3.2093 


1. 170 


G.()38 


1510 


1.921 


0, liGi 


,_i 


Hog 



Note. — [The complete version of this table is in the electronic edition of the Astrophysical Journal and from http://physics.uwyO.edu/agn/. The printed edition contains 
only a sample.] Units on are 10~ "'^^ erg cm~^ s~ "'^ A~ ■'^ . Units on are 10*^^ erg s""'^ . Both and \ are taken at 3000 A. Units on the Mg ii FWHM are km s~"^ . 

Units on A^gjj are 10^ ^0- 
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